Attractive and repulsive contributions of medium fluctuations to nuclear superfluidity 



o 
o 

(N 
tin 

cn 
<N 



> 

o 

<N 
O 

in 
o 



X 



G.GorP'^ F.RamponP'^ F.Barranco^ P.F. Bortignon^'^ R.A. Broglia^'^'^, G. Col6^'^ E.Vigezzi^ 
^ Dipartimento di Fisica, Universita degli Studi di Milano, via Celoria 16, 20133 Milano, Italy. 
^ INFN, Sezione di Milano, via Celoria 16, 20133 Milano, Italy. 
^ Departamento de Fisica Aplicada III, Escuela Superior de 
Ingenieros, Camino de los Descubrimientos s/n, 41092 Sevilla, Spain. 
^ The Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen 0, Denmark. 

(Dated: February 9, 2008) 

Oscillations of mainly surface character (5* = modes) give rise, in atomic nuclei, to an attractive 
(induced) pairing interaction, while spin (S = 1) modes of mainly volume character generate a 
repulsive interaction, the net effect being an attraction which accounts for a sizeable fraction of 
the experimental pairing gap. Suppressing the particle- vibration coupling mediated by the proton 
degrees of freedom, i.e., mimicking neutron matter, the total surface plus spin-induced pairing 
interaction becomes repulsive. 
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A central issue in a quantitative description of super- 
conductors (metals, doped fullerides, etc.) as well as su- 
perfiuid Fermi systems (^He, trapped gases of fermionic 
atoms, atomic nuclei, neutron stars, etc.) is related to 
the interaction acting between fermions and giving rise to 
Cooper pairs. The glue holding together the fermions of 
each Cooper pair, is the result of the bare interaction act- 
ing between fermions, strongly renormalized by medium 
polarization effects: Coulomb interaction plus plasmon 
and phonon exchange in metals 1], Van der Waals inter- 
action plus spin and density modes exchange in trapped 
gases of fermionic atoms 0, , strong force and spin and 
density modes exchange in the case of atomic nuclei and 
of neutron stars 0, H, Q • 

While broad consensus exists concerning the mecha- 
nism of electron-electron and electron-phonon interaction 
leading to superconductivity in metals, the situation is 
much less clear in the case of strongly interacting par- 
ticles in the different scenarios found in nature. In par- 
ticular, it has been found that medium polarization ef- 
fects associated with the exchange of vibrations lead to a 
quenching of the bare pairing interaction in the chan- 
nel in the inner crust of neutron stars 7] while, at the 
same time, account for a sizeable fraction of the pairing 
gap in open shell nuclei 8, 9, 10]. In the present paper 
we present, for the first time, evidence which allows to 
understand these seemingly contradictory results. 

We address the question at hand within the scenario 
provided by the paradigmatic (superfluid) open shell nu- 
cleus ^^^Sn. The starting point corresponds to the calcu- 
lation of the mean field potential and associated quasi- 
particle properties within the framework of Hartree-Fock 
plus BCS theory 11] using the SkM* force 12]. The 
polarization quanta were worked out within the frame- 
work of the quasi-particle random phase approximation 
(QRPA). The particle- hole residual interaction is derived 
in a self-consistent way from the Skyrme energy func- 
tional, with the exception of the spin-orbit and of the 
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FIG. 1: Generalized Landau parameters associated with the 
interaction SkM* defining the strength of the particle-hole in- 
teraction in the isoscalar (Fo), isovector (Fq), spin isoscalar 
(Co) and spin isovector (Gq) channels. In the inset the func- 
tions Fo + Fq (n-n interaction), Fo — Fq (n-p), Go + Gq (n-n) 
and Go — Gq (n-p) are also shown. 



Coulomb part (cf. 1^ for more details). On the other 
hand, we neglect the momentum-dependent part of the 
interaction in the calculation of the particle-vibration 
coupling discussed below 14]. The relevant part of the 
the particle-hole interaction can then be written 



Vph{r,r) = S{r- r) x 
x{[i^o + i^or-r'] + [(Go 



■G',T-T')a-a']}. (1) 



We shall only consider the -r^ term, in keeping with the 
fact that we are here interested in the neutron-neutron 
pairing interaction. Off-diagonal terms are associated 
with charge-exchange modes. Thus, in lowest order, they 
do not contribute to the neutron-neutron interaction, but 
are expected to be of relevance in the discussion of the 
proton-neutron pairing interaction. 
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The functions Fo{r),F^(rY Go{r) and G[)(r) (gen- 
eralized Landau-Migdal[l3 [l^ parameters) controlling 
the isoscalar and isovector (spin-independent and spin- 
dependent) channels are displayed in Fig. Q 

Strictly speaking, in the case of atomic nuclei spin 
is not a good quantum number with which to identify 
the polarization quanta, because of the strong spin-orbit 
term present in these systems. We have thus adopted the 
criterion of distinguishing between natural (tt = (— l)'^) 
and non-natural (tt = —(—1)'^) parity modes, where J in- 
dicates the total angular momentum of the quanta. Vi- 
brations of multipolarity and parity = 1+, 2+,2~, 
3+,3~, 4+,4~,5+ and 5~ were worked out. Those hav- 
ing energy <30 MeV were used in the calculation of the 
induced interaction (cf. Fig. Efa)). From the QRPA cal- 
culation one gets , together with the energy of the ex- 
cited states, their transition densities which will be used 
as form factor for the particle- vibration coupling vertex 
(cf. Fig. Hb)): 





/ 










1 


\ 




(a) 


V 


Phonon with, natural parity 


^ Vph 










Phonon with non-natural parity 


/ " '"I 












(b) 





727+1 



FIG. 2: (a) diagram depicting the pairing interaction induced 
by the exchange of phonons; (b) particle-vibration coupling 
vertex making it explicit the dominant part of the particle- 
hole interaction giving rise, through the sum of bubble dia- 
grams, to the corresponding QRPA modes (wavy line). 
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similar way, the vertices produced by the spin dependent 
part of the residual interaction are 
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,V2^^-) J )) 



■Uu2Vui) X {ui\\i^[YL X cr]j\\u2)(Pui{r)(Pu2{r). 

(3) 

The function Q is associated with the response to ex- 
ternal fields which induce a density oscillation; it vanishes 
for phonons of unnatural parity. The function is in- 
stead associated with the response to magnetic external 
fields and hence the coupling to excited states mediated 
by the part of the residual interaction which depends on 
the spin, and it applies to phonons of both unnatural 
(when J L) and natural (when J = L) parity. The 
index i labels the different vibrational modes of a given 
spin and parity, while X and Y are the forwardsgoing and 
backwardsgoing QRPA amplitudes of the corresponding 
modes, and the index u denotes the quantum numbers 
n, j of the single particle states. 

By using Eq. (^P), one can find the expression of the 
vertices produced by the spin independent part of the 
residual interaction: 



u m 



{v'm'l [Fo(r) + F^(r)f • f] S{r- r')\iym; J" Mi). (4) 

It can be rewritten by using the multipole expansion for 
the 6 function as 
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^p^j-K^ and Spljnp being the neutron and proton contribu- 
tions to the transition densities defined in Eq. Q. In a 
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{u'm'l [Go{r) + G[)(r)r • r V ' ^'S{r- i^)\vm; r 
and, as before, they can be expanded in the form 
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dnp., [(Go + G'o)5p^.^„ + (Go - G'o)5pij.r.p]M7) 



where ^Pjtt^^ and Spj-j^i^p are respectively the neutron 
and proton contributions to the transition densities de- 
fined in Eq. 

These particle-vibration coupling matrix elements, to- 
gether with the energies of the QRPA modes and the HF 
single-particle energies, are the basic ingredients needed 
to calculate the pairing induced interaction Vjnd (cf. Fig. 
I2fa)) within the framework presented in ref. UQj. The 
matrix elements between two pairs of neutrons in time- 
reversal states are given by 

< v' m v' m'\vind\^'^^'^ >= 

(5j 



Eq — E, 



'int 



where the sum is over the phonons of multipolarity , M 
obtained in the QRPA calculation. In the denominator, 
Eint = {\ey-eF\^\ey' -eF\^huOi) denotes the energy of 
the intermediate state given by two particles (whose en- 
ergy is calculated respect to the Fermi energy ep) and one 
vibration, while Eq is the energy of the correlated two- 
particle state, which must be obtained self-consistently: 
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FIG. 3: Diagonal induced pairing matrix elements resulting 
from the exchange of phonons with natural parity (filled cir- 
cles) and those resulting from the exchange of phonons with 
non-natural parity vibrations (empty circles), displayed as a 
function of the energy of the single-particle state . 



the denominator is therefore always negative. The impor- 
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FIG. 4: (a) The state dependent pairing gap as a function of 
the single-particle energies obtained by solving the BCS equa- 
tions associated with the total induced interaction matrix el- 
ements; (b) same as (a) but for the matrix elements produced 
including only the spin independent part of the particle-hole 
interaction. 



tant conclusion is that the sign of the matrix element (O 
then depends on the relative magnitudes of the attrac- 
tive contribution arising from the spin- independent 
term and of the repulsive contribution —g'^ arising from 
the spin-dependent term. 

The matrix elements Q are then used to obtain the 
ISn neutron pairing gap, solving the BCS gap equation 



A,. = 



where 
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(9) 



V(2j + l)(2j' + 1) < jj; 0|vi„d|i'i'; > (10) 

The two particle wavefunction |j'/;0 >, coupled to 
zero angular momentum, contains an admixture of sin- 
glet and triplet components, with about equal weight. It 
is well known (see e.g. IIS l2l| ) that in infinite matter, 
when only a species is present, like in neutron matter, the 
matrix elements associated with the singlet component 
have a simple character, depending on the spin S of the 
exchanged fluctuation: for density modes, characterized 
by S =0, the matrix elements are negative, while for spin 
modes (5* = 1) they are positive. The suppression of the 
pairing gap produced by the bare neutron-neutron force 
through medium polarization effects in neutron matter is 
associated with the dominance of the repulsive contribu- 
tion from the spin modes over the attractive contribution 
from the density modes. 

In the case of ^^^Sn, we observe that the diagonal 
matrix elements < jj]^\vind\jj]^ shown in in Fig. 
13 are systematically attractive or repulsive, depending 
on whether the exchanged fluctuations are respectively 
of natural or unnatural parity, in direct correspondence 
with the case of infinite matter. On one hand, this is 
because for non natural parity modes only the spin de- 
pendent vertices Q, which have a 5* = 1 character, con- 
tribute, while for natural parity modes, the spin indepen- 
dent matrix elements are the dominant ones, and one 
can show that they are the only ones contributing to the 
diagonal matrix elements. On the other hand, the con- 
tributions of the matrix elements from triplet component 
of the two-particle | jjO > wavefunction are small and of 
variable sign [2^ . 

However, in constrast with neutron matter, the result- 
ing total matrix elements are predominantly attractive 
(in any case around the Fermi energy) 23]. 

The resulting state dependent pairing gap obtained by 
solving the BCS gap and number equations making use of 
the (total) induced pairing matrix elements is depicted in 
Fig. Efa). For states close to the Fermi energy the gap ac- 
counts for a consistent fraction of the experimental value 
(1.4 MeV). If one solves the BCS equations considering 
only the exchange of density modes (i.e. neglecting the 
contributions from Eq. Q), one obtains values which 
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FIG. 5: (a): the dashed and dot-dashed curves are respec- 
tively the neutron and proton transition densities associated 
with the 2^ phonon with energy 1.5 MeV while the solid curve 
is the wave function of the lhii/2 state (in arbitrary units), 
(b): same as (a) but for the 3^ phonon with energy 4.35 MeV 
and the 2d3/2 state. 



are, in average, larger (cf. Fig. Efb)). In fact, the ex- 
change of = 1 modes quenches the pairing gap arising 
from the exchange of only 5 = modes by roughly 30% 

To gain insight about the peculiar features of finite 
nuclei, as compared to the case of infinite systems, it 
is useful to study the radial dependence of the particle- 
vibration coupling vertices shown in Fig. El^b). The in- 
duced pairing matrix elements associated with natural 
parity modes have a clear surface character (cf. ref 25]). 
In particular this is the case for the most attractive pair- 
ing matrix element which is associated with the ^^xxj^^^ 
(^1/^11/2 = -8-0'' MeV, eF = -8.50 MeV) configuration 
(cf. Fig. Because of its large centrifugal barrier, 

the wave function of this single-particle state is mainly 
concentrated at the nuclear surface. The main contribu- 
tion to the corresponding induced pairing matrix element 
arises from the exchange of a 2+ phonon (of energy 1.5 
MeV) between the two nucleons moving in time reversal 
states in the hii/2 orbital. The associated proton and 
neutron transition densities depicted in Fig. ISI^a) tes- 
tify to the fact that this phonon has the character of a 
surface vibration. Concerning the most repulsive matrix 



elements, we have found that the corresponding unnatu- 
ral parity phonons are volume modes. In particular, one 
of the largest (positive) matrix element is associated with 
the 2d3^2(0) configuration (62^3/2 =-8-52 MeV). Because 
of the low angular momentum, one finds that a consistent 
fraction of the corresponding wave function is concen- 
trated in the interior of the nucleus. This state can thus 
couple efficiently with phonons of volume character. In 
fact, the major contribution to the corresponding matrix 
element is due to the exchange of the 3+ vibration (with 
energy at 4.35 MeV) which is a mode with a large volume 
component as testified by the corresponding proton and 
neutron transition densities shown in Fig. IS^b). One can 
conclude that states lying close to the Fermi energy with 
high 2 and thus localized at the surface mainly feel the 
(attractive) coupling arising from the exchange of surface 
vibrations. The situation is expected to be quite different 
in the case of infinite neutron matter. In fact, in going 
from the finite to the infinite system the collectivity of 
the natural parity modes, mostly surface modes, will be 
strongly reduced, while not much is expected to happen 
to the volume modes. 

Furthermore, in going from nuclear to neutron mat- 
ter, many attractive contributions vanish. In fact, if 
we turn off the neutron-proton particle-hole interaction 
contributing to the basic vertices displayed in Fig. ^ a 
strongly net repulsive induced interaction is obtained (cf. 
Fig. a situation which much resembles the neutron 
star case. This result can be understood by realizing that, 
quite generally, the dominant contribution to the spin in- 
dependent (and therefore attractive) induced matrix el- 
ements arises precisely from the neutron-proton part of 
the particle-hole interaction, which is proportional to to 
(Fo - Fq)5pp (cf. Eq. (4) and Fig. QJ. The remaining 
part of the spin independent interaction depends on the 
function Fq + Fq (corresponding to the particle-phonon 
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FIG. 6: The diagonal matrix elements produced by the ex- 
change of phonons with natural parity (filled circles) and those 
produced by the exchange of phonons with non-natural par- 
ity (empty circles) when the proton part of the phonon wave 
function is not included in the calculation, are displayed as a 
function of the energy of the single-particle states. 
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coupling mediated by Spj-r,^), which is rather weak, and 
for the SkM* interaction adopted here, even displays a 
node at the nuclear surface. The induced interaction is 
then dominated by the (repulsive) spin-dependent matrix 
elements proportional to Go + Gq (corresponding to the 
neutron-neutron particle-phonon coupling mediated by 
^Pj-^Ln^ cf. Eq.(5)), which are large and without nodes. 

We have chosen to adopt a force like SkM*, which has 
been tested by different groups over the years. We are 
well aware that our results could change adopting a dif- 
ferent effective force, especially in connection with the 
fact that it is difficult to fix the spin-dependent part of 
the interaction, also because of the scarcity of experimen- 
tal constraints. In fact, different approaches have been 
adopted in the literature. In particular, 1+ modes have 
been calculated in the framework of the extended the- 
ory of finite Fermi systems, using constant values of the 
Landau parameters Go and Gq ^y, [23 • Recently, a new 
specific parametrization of the Skyrme type has been in- 
troduced, SkO', which takes into account the effects of 
time-odd spin-isospin couplings, and has been adopted 
for the description of the spin-flip transitions like the 
Gamow- Teller resonance |2^|23. This force has not been 
yet extensively checked in the non-charge-exchange chan- 
nel, and the large and negative value of Go associated 
with this force can lead to too strongly collective or un- 



stable solutions in the spin isoscalar channel. Even using 
different effective forces, however, the main qualitative 
aspects of our results are likely to remain true, as they 
are based on quite general features - like the dominance 
of neutron-proton over neutron-neutron interaction, and 
the stronger collectivity and surface localization of the 
low-lying density modes with respect to the spin modes. 

We conclude that the exchange of low-lying surface 
vibrations (in which neutrons and protons participate 
on equal footing) between pairs of nucleons moving in 
time reversal states close to the Fermi energy, leads to a 
sizeable attractive pairing interaction which accounts for 
about 70% of the pairing gap. The inclusion of spin (vol- 
ume) modes, reduces this contribution by 30% in the case 
of finite nuclei, bringing the induced pairing contribution 
to the pairing gap to a value of the order of ~50%, the 
other half coming from the bare nucleon-nucleon force. 
The attractive character of the effective interaction in fi- 
nite nuclei is found to be associated with the efficient 
coupling of the single-particle states lying close to the 
Fermi energy to collective surface vibrations, as well as 
with the contribution of the proton-neutron part of the 
particle- vibration coupling. Without these two elements, 
spin modes dominate and the effective interaction be- 
comes repulsive, as in neutron matter. 
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